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RESUMEN 

Se revisan trabajos sobre la morfologia y distribucion galactica de las nebulosas planetarias (PNs), asf como 
los avances recientes en modelaje MHD de las PNs. Llegamos a una explicacion tentativa sobre la conexion 
entre las clases morfologicas y su distribucion galactica. 

ABSTRACT 

We review recent studies on the morphology and distribution with respect to the galactic plane of planetary 
nebulae (PNs), as well as recent advances in MHD modeling of PNs. We discuss a tentative explanation for 
the connection between morphological classes and their galactic distribution. 
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1. INTRODUCTION 

The connection between PN morphologies and 
their distribution with respect to the midplane of 
the Galaxy has been studied, with notorious ad- 
vances, during the last three decades. The search 
for systematic segregations among PNs of different 
shapes started with the analysis by Greig (1972). 
He had two main morphological classes, binebulous 
(50%) and circular (50%), and found that the bineb- 
ulous class had a higher strength in the forbidden 
lines [O HI], [O II] and [N II]. He also found a hint 
of a correlation between this class and smaller dis- 
tances to the galactic plane on average. Later, us- 
ing the abundance of chemical elements, Peimbert 
(P), Torres-Peimbert (T-P) and collaborators (PT- 
P71, T-PP77, P78, T-PP79, PT-P83, T-PP83; see 
T-PP97 for a review) made a PN classification with 
four classes, from Type I to Type IV, with decreas- 
ing abundances of Helium and heavy elements. They 
found that most Type I PNs have bipolar shape, and 
this lead to the conclusion that there might be pop- 
ulation differences between bipolar and circular PNs 
in our Galaxy. In the abstract of their seminal pa- 
per, Calvet and Peimbert (1983) wrote: "/t is sug- 
gested that the bipolar nature of PN of Type I can be 
explained in terms of their relatively massive progen- 
itors (M > 2.4 Mq ), that had to lose an appreciable 
fraction of their mass and angular momentum dur- 
ing their planetary nebula stage. In a first mass-loss 
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stage at low velocity it is produced a disk and in a 
second mass-loss stage at higher velocity the matter 
is limited by the disk, giving rise to the bipolar struc- 
ture" . 

Later, several authors made improvements in our 
knowledge of this issue. Zukerman & AUer (1986) 
and Zukerman & Gatley (1988) classified 108 PNs in 
bipolar, round, disk-like, and annular. They found 
an anti-correlation between metallicity and Galactic 
latitude, as well as a segregation of morphological 
types according with Galactic latitude, in agreement 
with previous studies. 

Balick (1987) made a major contribution to the 
morphological classification (50 PNs), and named 
the classes as round, elliptical and butterfly. At the 
same time, Chu et al.(1987) (126 PNs) pubhshed a 
catalog of PNs with more than one shell (50 % of the 
sample), and named this new class as multiple-shell 
PNs. 

During the last decade, two important surveys 
were carried out on both hemispheres. For the south- 
ern hemisphere. The ESO Survey, the study was 
published in a series of papers by Schwarz, Corradi 
& Melnick (1992), Stanghellini, Corradi & Schwarz 
(1993) (250 PNs), Corradi k Schwarz (1995), and 
Corradi (2000) (400 PNs). For the northern hemi- 
sphere. The I AC Survey, the study was published by 
Manchado et al.(1996) (243 PNs) and Manchado et 
al.(2000) (255 PNs). The first detailed study of the 
differences between ellipticals and bipolars was done 
by Corradi & Schwarz (1995) with the ESO Survey 
(see also Corradi 2000). They found that the bipo- 
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lar class has a smaller scale height, 130 pc, than the 
one for ellipticals, 320 pc. Also, bipolars have the 
hottest central stars among PNe, and display smaller 
deviations from pure circular Galactic rotation than 
other morphological types. In addition, bipolars also 
display the largest physical dimensions and have ex- 
pansion velocities of up to an order of magnitude 
above the typical values for PNe. These properties, 
together with the chemical abundance results by Cal- 
vct & Peimbert (1983), indicate that bipolar PNe are 
produced by more massive progenitors than the re- 
maining morphological classes. 

The distributions inferred by these surveys rep- 
resent the main motivation for the present paper, 
which discusses these results in the light of the new 
MHD models. 

2. BIPOLARS: ROTATION AND ANGULAR 
MOMENTUM 

If bipolarity truly comes from the stellar rotation 
and the loss of angular momentum in massive PN 
progenitors, one can learn from the case of the most 
massive stars and their associated nebulae. In the 
Hertzsprung-Russell diagram, LBV stars like rj Car 
lie close to the upper limit in temperature and lu- 
minosity beyond which no normal stars are observed 
(the Humphreys-Davidson Limit), near the location 
of the Eddington limit. 

A massive star expanding at a constant luminos- 
ity L would — in the non-rotating case — reach the 
Eddington limit due to an opacity increase in the 
surface layers. The Eddington limit is reached when 
F = L/L-Edd = 1) where the Eddington luminosity 
is I/Edd = 4:TrcGM/K, with c the speed of light, M 
and L the mass and luminosity of the star, and k the 
opacity at the stellar surface. Note that when F — > 1, 
the escape velocity v^sc = [(1 - F)2GM/i?]V2 ^ o. 
This implies a very slow wind velocity (~ 0) at the 
Eddington limit. Thus, when a star passes through 
a close approach to the Eddington limit, its stellar 
wind will follow a fast sloiv fast sequence. If we 
now include rotation, the critical rotational veloc- 
ity also follows the same behavior. The combina- 
tion of radiation pressure and centrifugal force ex- 
ceeds surface gravity at the equator when iicrit = 
[(1 - F)GM/i?] and v^it ^ when F ^ 1. How- 
ever, given that the rotational velocity v^ot can be 
small but is always different from zero, critical ro- 
tation occurs before the Eddington limit is reached. 
This is the so called limit (Langer 1997). Winds 
with mass loss rates large enough to halt and re- 
verse the stellar expansion are thought to occur when 
^ = ^^rot/^^crit approaches unity (Friend & Abbott 



1986). Thus, rotationally driven outflows will always 
occur before the star reaches its Eddington limit. 

Using this idea, Langer et al. (1999a) simu- 
lated the LBV outburst phenomenon by allowing the 
star to pass through three successive stages — pre- 
out burst, outburst and post-outburst — with cor- 
respondingly different stellar winds. To include the 
effects of stellar rotation on the winds, they used 
the analytic model of Bjorkman & Cassinelli (1993), 
with stellar parameters appropriate to r] Car. The 
models were not intended to exactly fit all proper- 
ties of the Homunculus nebula. It is however strik- 
ing how well these models, without much fine tun- 
ing, not only reproduce the large-scale, bipolar mor- 
phology, but also the small scale turbulent structure 
seen in high-resolution observations of the Homuncu- 
lus. This empirical result, together with the observa- 
tional facts that all LBV nebulae are bipolar (Nota 
& Clampin 1997), strongly indicates a link between 
rotation and bipolarity. 

What about Bipolar PNs? In order to evaluate 
whether or not AGB rotation rates can produce situ- 
ations similar to the one discussed above, one needs 
to derive F and the rate of critical rotation. Di- 
rect determinations of rotation rates for AGB stars 
are rare and mostly provide upper limits, of the or- 
der of a few km s^^ (except for V Hydrae with 
vsini ~ 13 km ). Using the available informa- 
tion for main sequence stars in the initial mass range 
associated with PN progenitors (from ^ 0.8 Mq to 
about 5Mq) these can be divided into two groups. 
Stars with masses below ^ 1.3 M0, with subpho- 
tospheric convective zones, spin down during their 
main sequence evolution due to flares and magne- 
tized winds. When they reach the AGB phase, their 
resulting rotation velocities are below 0.01 km s^^ . 
Thus, the expected rotation speeds of PN progen- 
itors with masses below ~ 1.3 M©) should be very 
small. 

Stars above ~ 1.3 Mq do not have convective en- 
velopes during core hydrogen burning, and they ap- 
pear to remain as rapid rotators throughout their 
main sequence evolution. When they develop sur- 
face convection on their way to the red giant branch, 
they create a hydrogen burning shell which separates 
the helium core from the envelope. The helium core 
evolves decoupled from the envelope and retains its 
angular momentum, i.e. the entropy barrier of a nu- 
clear burning shell prevents that angular momentum 
can leak out of the core. Also in this case, how- 
ever, the angular momentum of the hydrogen enve- 
lope will be lost, either due to magnetic braking or 
to mass loss and reexpansion of the convective enve- 
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lope on the AGB (cf. Langer et al. 1999b). As these 
stars move to the thermally pulsing AGB phase, 
the H and He burning shells periodically switch on 
and off. Thus, the presumable barriers vanish peri- 
odically and core-envelope angular momentum ex- 
change can occur during this stage. Such an ex- 
change is very likely, since mixing of matter through 
the core boundary is known to occur in this phase, 
which is necessary to activate the ^^C{a,n) reac- 
tion, which is the neutron source to operate the s- 
process. To estimate the resulting rotational velocity 
of the envelope (of about 0.1 Mq) at the end of the 
AGB phase, we may approximate the specific angu- 
lar momentum of the core by its main sequence value, 
which is transferred to the envelope. For a core 
of ^ 0.5 and a main sequence radius of about 
O.IR© we get j ~ 10^^...10^* cm^ s~^ Adopting an 
average radius of the AGB envelope of lOOR©, we 
obtain a rotation velocity of Vrot — 10^^'^ cm^ s~^ x 
0.5Mq/100Rq/0.1Mq ~ 2 km s'^ . Regardless of 
how crude this estimate may be, the net result is that 
stars at the tip of the RGB or AGB phases may be 
subject to a significant spin-up of their surface lay- 
ers (Heger & Langer 1998). Therefore, we propose 
here that values of = VroL/vcrU close to 1 may 
be appropriate for AGB single stars above ~ 1.3 Mq 
during the phase of PN ejection. We parametrize the 
"superwind" in the same way as the giant LBV out- 
bursts (Eddington parameter F = L/Ledd close to 1; 
see Garcia-Segura et al. 1999 for further details). 

There arc variations of this scenario and the pos- 
sibility exists that stars which are not able to reach 
the f2 limit on their own might be able to do so in 
the presence of a companion. For example, main 
sequence wide binary systems, which become close 
binaries (either detached or attached) at the AGB 
phase, can be subjcc;t to a very effective; spin-up by 
their companion. Depending on the binary separa- 
tion and masses, the spin-up can be produced by 
tidal forces (slow process) or by spiral-in (relatively 
fast process). In either cases, the orbital angular 
momentum is transfered to the star. 

Sokcr (1995) has studied the tidal spin-up ap- 
plied to the formation of elliptical, but the theory can 
be also applied to the formation of bipolars (Soker 
& Rappaport 2000). Extreme cases of the same 
scenario is the common envelope evolution scenario 
(see for example Livio & Pringle 1996 and references 
there in, also Reyes- Ruiz & Lopez 1998). 

Thus, a binary system of such a class can form a 
bipolar nebula, however, the nebula will show point- 
symmetric features, as we will discuss later on sec- 
tions 4 and 5. Also, these nebulae do not have to 



show an enrichment by heavy elements, and so, they 
do not necessary fit into the Type I Peimbert cate- 
gory. Finally, since their primary stars may be below 
^ 1.3 , their scale height over the plane can be 
much larger than the one for classical Type I Peim- 
bert nebulae (see section 5). 

3. ELLIPTICALS: A MHD PHENOMENON 

Axisymmetric flows can be produced by a mag- 
netized wind with or without the existence of equa- 
torial density enhancements (EDE) (see Rozyczka & 
Franco 1996, Garcia-Segura 1997, Garcia-Segura et 
al. 1999 and Garcfa-Scgura & Lopez 2000). An EDE 
can be formed with a small amount of rotation at the 
AGB phase for example (Ignace et al.l998) or by a 
dipolar magnetic field (Matt et al.2000). The mag- 
netic field at the surface of a post- AGB star can be 
transported out by its wind, similar but not equal to 
the solar case. Because of stellar rotation, the mag- 
netic field in the wind is dominated by a toroidal 
component. The resulting toroidal field has a mag- 
netic tension associated with it. Thus, the general 
effect of the magnetic tension is the elongation of 
the nebula in the polar direction. The mechanism re- 
sponsible for the elongation is described by Rozyczka 
& Franco (1996) and Franco et al. (2001, this con- 
ference) in detail, and successful examples of the for- 
mation of jets and ansae (FLIERS) can be seen in 
all above articles at the beginning of this section. 

Unfortunately, there is no information about stel- 
lar magnetic fields at the surfaces of post- AGB stars, 
which is in one of the most interesting regions of the 
HR diagram, between the AGB phase and the white 
dwarf phase. Rcid et al.(1979) found strong circular 
polarization in the OH masers from U Ori and IRC 
-1-10420. The OH masers indicate fields of about 10 
milligauss emanating from regions of about 10^''' cm 
from the star. They suggest that the field strengths 
at the stellar surfaces are on the order of 10 and 
100 Gauss, respectively. White dwarfs arc known to 
have fields of up to 10^ - 10^ Gauss (Schmidt 1989). 
So, stellar field intensities between the wide range 
defined by the solar value {B^ — 2 Gauss, in av- 
erage) and 10* Gauss are acceptable for modeling. 
Note that observations for the detection of magnetic 
fields in the swept- up shells of PNc (eg Terzian 1989) 
have been tried, but however, they give us very little 
information about the field intensity of the nucleus, 
since that diluted field (if detected) will correspond 
to the already expanded and compressed AGB gas. 
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4. THE POINT-SYMMETRIC SUB-CLASS: 
SIGNATURES OF BINARIES 

A particularly intriguing case in PNs morpholo- 
gies are those that display point-symmetric struc- 
tures. At first view, the point-symmetric morpho- 
logical class does not look very important. But, a 
careful inspection to the statistically, complete sam- 
ple of the lAC morphological catalog (Manchado et 
al.l996) reveals 40 objects with some degree of point- 
symmetric features, which represent 19% of the total 
list (215) with a well defined morphology. Note that 
many of these nebulae are not classified as point- 
symmetric. In fact, they have been well classified 
inside the categories of bipolars and ellipticals, such 
as the case of the Dumbel nebula. This fact points in 
the direction that point-symmetry is a common fea- 
ture related to any morphological class, instead of a 
separate group (Manchado et al.2000; see also Guer- 
rero, Vazquez & Lopez 1998). Such a large fraction 
of the sample (19%) suggest that the reason which 
produces point-symmetry should be very common 
indeed. 

The most convincing solutions up to now, for 
the formation of point-symmetric nebulae require 
the existence of a binary system, and a magneto- 
hydrodynamical coUimation of the wind, cither for 
accretion disk winds (Livio & Pringle 1996; see also 
Mastrodemos & Morris 1998, Reyes-Ruiz & Lopez 
1998 and Blackman et al. 2001) or for stellar winds 
(Garcia-Segura 1997, Garcia-Segura & Lopez 2000). 

Since the toroidal magnetic field carried out by 
the wind, cither stellar or coming from a disk (for 
the last one see Contopoulos 1995), is always perpen- 
dicular to the rotation axis of the central star/disk 
(see Figure 3 in Garcia-Segura ct al.2GG0), any kind 
of misalignment from the axis (wobbling instabil- 
ity, precession, steady tilt respect to the equato- 
rial density enhancement) will produce "naturally" 
a point-symmetric nebula. In either case, it is easy 
to imagine the topology of the magnetic field lines 
in such scenarios, i.e., multiple rings centered and 
aligned along the "time-dependent" spin axis of the 
star/disk. 

Close binaries arc a necessary condition in Livio 
& Pringle (1996) , while wide binaries are a suffi- 
cient one in Garcia-Segura (1997) and Garcfa-Segura 
& Lopez (2000). As described in section 2, close bi- 
naries, either attached or detached, will be in favor 
of forming bipolars with point-symmetric features, 
while wide binaries will be in favor of ellipticals with 
point-symmetric features. 



5. DISCUSSION 

The average scale height over the plane, < z >, 
for different PN morphological classes can be com- 
pared with those of different stellar masses and pop- 
ulations. For instance, we know that massive stars 
are located much closer to the galactic plane than 
the population of stars with lower initial mass. The 
outcome from the ESO and lAC surveys is certainly 
coincident: the bipolar (B) class has < z > = 130 pc 
(ESO) and < ^ > = 179 pc (lAC), for ellipticals (E) 
< z > = 320 pc (ESO) and < z > = 308 pc (lAC), 
and for rounds (R) < z > = 753 pc (lAC). 

In the recent analysis of the lAC Survey by Man- 
chado et al. (2000), the bipolar (BPS) and eUiptical 
(EPS) objects with point-symmetric features were 
separated from those which do not present such a 
kind of symmetries, i.e., from the B and E classes 
respectively. The new results from the lAC Survey 
are given in the Table 1. 

Comparing the results of the lAC Survey with 
those described by Miller & Scalo (1979) for the aver- 
age scale heigth of stars with different masses, an av- 
erage value equal or smaller than 110 pc corresponds 
to stars with initial masses above 1.9 M© . Lower 
mass stars have average scale heiths well above this 
value. These results arc in line with the pioneer- 
ing suggestion of Calvet & Peimbert (1983) and the 
more recent discussion made by Garcfa-Segura et 
al.(1999). 

To summarize, from the new MHD studies (sec- 
tions 2, 3 and 4) we can conclude that: 

Bipolarity 4^ Omega Limit Angular Momentum 
(either orbital or stellar) 

EUipticity <^ MHD effects ^ Symmetric Mi- 
crostructures (FLIERS, jets, ...etc) 

Point- Symmetry Binaries Precession, Wob- 
bling, ... 

TABLE 1 



Morphological Class < z > 

1«_ 

B 110 
BPS 248 
E 308 
EPS 310 
R 753 
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< Z > The lAC Catalog 




Fig. 1. 

Finally, the relation between morphology and 
galactic distribution could be explained as (see Fig- 
ure 1): 

• Bipolars Type I Peimbert = (B): 

Small < z > (110 pc) <s=^ Massive Progenitor 



■^=> Stellar Rotation • 
Bipolarity 



O Limit 



Classical 



• Tidal Bipolars = (BPS): 

Moderate < z > (248 pc) <^=^ Non-Massive 
Progenitor in Close Binary System <J=> Tidal 
Spin-Up <S=^ Shaping by fl Limit + MHD 
Effects -|- Precession/ Wobbling <s=^ Bipolarity 
with Point-Symmetry (Lobes, FLIERS, Jets) 

• Ellipticals = (E): 

Medium < z > (308 pc) Non-Massive 
Progenitor <J=^ Shaping by MHD Effects ■^=^ 
FLIERS & Jets with Axisymmetry 

• Ellipticals with Point-Symmetry = (EPS): 

Medium < z > (310 pc) <^=^ Non-Massive Pro- 
genitor in Wide Binary System <==^ Shaping by 
MHD Effects + Precession FLIERS & Jets 
with Point-Symmetry 

• Round = (R): 

Large < z > (753 pc) «^=> Low-mass Progenitor 
<s=^ Neither Rotation nor MHD effects 
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